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We  review  several  experimental  methodologies  for  the  determination  of  nonlinear  absorption 
coefficients  of  organic  chromophores  in  either  solution  or  thin  film  form.  The  goal  here  is  to  also 
be  able  to  identify  and  quantify  the  nonlinear  absorption  process  or  processes,  e.g.  the  essentially 
instantaneous  two-photon  absorption  (2PA)  or  the  cumulative  excited-state  absorption  (ESA). 
These  methods  include  the  direct  methods  of  nonlinear  transmission  and  Z-Scan,  along  with  the 
indirect  methods  of  two-photon-induced  fluorescence,  excite-probe  approaches,  and  Degenerate 
Four-Wave  Mixing.  We  also  discuss  both  frequency  degenerate  and  frequency  nondegenerate 
processes  with  emphasis  on  2PA. 

It  is  important  to  note  that  a  single  method  is  generally  incapable  of  both  measuring  and  identifying 
the  source  of  the  nonlinearity,  thus  it  is  usually  necessary  to  employ  more  than  a  single  method  to 
determine  the  sign,  magnitude  and  physical  processes  involved.  For  example,  single  beam 
experiments  cannot  give  the  temporal  response  of  the  nonlinearity  while  excite-probe  methods  do; 
however,  single  beam  measurements  lend  themselves  better  to  absolute  calibration.  In  addition 
we  include  information  on  nonlinear  refraction,  which  always  accompanies  nonlinear  absorption 
processes.  For  example,  the  nondegenerate  bound-electronic  nonlinear  refraction  and  ultrafast 
nondegenerate  nonlinear  absorption  come  from  the  real  and  imaginary  parts  of  the  nondegenerate 
third-order  nonlinear  susceptibility  and  are  related  by  Kramers-Kronig  relations  via  causality. 
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1.  INTRODUCTION: 

In  this  chapter,  we  describe  several  techniques  for  measuring  nonlinear  absorption  (NLA)  with 
special  emphasis  on  the  instantaneous  nonlinear  response  of  two-photon  absorption  (2PA)[1]  due 
to  its  extensive  use  in  microscopic  imaging  and  3D  stereolithography/additive  manufacturing.  [2] 
We  also  give  some  information  on  nonlinear  refraction  (NLR)  measurements  since  NLR  almost 
always  accompanies  NLA  signals.  In  general,  multiple  experimental  techniques  are  needed  to 
fully  unravel  the  various  physical  processes  involved  in  the  nonlinear  absorption  for  a  given 
material,  as  it  is  rare  to  find  only  a  single  dominant  nonlinear  response.  For  example,  2PA  not  only 
causes  a  change  in  transmission  but  also  results  in  the  generation  of  excited  states.  These  excited 
states  can  also  change  the  absorption  spectrum  of  the  material  which  can  easily  be  confused  with 
the  2PA  itself. 

Additionally,  NLR  accompanies  NLA  and  transforms  into  a  spatial  redistribution  of  energy  upon 
propagation.  Thus,  in  any  experiment  to  monitor  transmittance  changes,  considerable  caution 
needs  to  be  taken  to  ensure  that  the  entire  beam  is  collected  by  the  detector  for  all  input  irradiances. 
Otherwise  the  NLR  may  be  easily  mistaken  for  NLA.  Also,  in  order  to  be  able  to  separate  NLA 
and  NLR,  the  conditions  for  external  self-action,  i.e.,  thin-sample  approximation,  need  to  be 
satisfied  putting  stringent  requirements  on  the  experimental  parameters. [3]  This  will  be  discussed 
in  more  detail  later.  Even  when  NLA  and  NLR  can  be  separated,  obtaining  reliable  values  for 
nonlinear  coefficients  can  be  difficult.  For  example,  if  the  NLA  is  irradiance  dependent  (as  for 
ultrafast  nonlinearities)  errors  come  from  energy  calibration  (e.g.  +5%),  temporal  measurements 
(e.g.  ±10%)  and  spatial  energy  distribution  measurements  (e.g.  ±10%  for  Gaussian  beams  but  for 
two  dimensions).  These  errors  add  along  with  fitting  errors  (e.g.  ±10%)  to  give  an  overall 
confidence  error  for  the  third-order  nonlinear  response.  The  example  errors  given  in  the  previous 
sentence  add  to  give  a  total  error  for  the  NLA  or  NLR  coefficient  of 

±Vo.052  +0.12  +2x0.12  +0.12  =±21%  .  Measurement  of  higher-order  responses  would  result  in 
larger  errors.  All  of  these  difficulties,  as  well  as  light  simply  missing  the  detector  and  incorrect 
identification  of  the  physical  processes,  have  led  to  reported  values  of  nonlinear  coefficients  that 
vary  over  orders  of  magnitude  from  one  publication  to  the  next.  [4-7]  While  such  large 
discrepancies  are  becoming  rarer  due  in  part  to  more  reliable,  stable  and  well  characterized  optical 
sources  and  to  our  more  complete  understanding  of  the  nonlinearities  themselves  along  with  the 
potential  measurement  problems,  obtaining  reliable  coefficients  with  small  confidence  errors  is 
still  challenging. 

2.  DESCRIPTION  OF  NONLINEAR  ABSORPTION  AND  REFRACTION  PROCESSES: 
2.1  2-Photon  Absorption  and  bound-electronic  nonlinear  refraction 

The  second  term  after  the  linear  optics  term  in  Bloembergen’s  [8]  expansion  of  the  electric 
polarization 

P  ~  £q(x{1)E  +  X(2)E2  +  T(3)£3  +...),  (1) 

is  the  third-order  polarization,  which  is  characterized  by  the  third-order  nonlinear  susceptibility 
y3).  This  is  the  lowest-order  term  in  the  expansion  that  leads  to  “self-action”  nonlinearities,  i.e. 
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nonlinear  absorption,  which  is  proportional  to  the  imaginary  part  of  ^3>  and  nonlinear  refraction 
proportional  to  the  real  part  of  y3). 

Assuming  that  the  sample  is  thin  in  that  the  sample  thickness  is  less  than  the  Rayleigh  length  and 
the  nonlinear  phase  does  not  propagate  far  enough  to  cause  an  irradiance  change  within  the  sample 
(i.e.  external  self-action[3]),  the  resulting  equations  for  the  propagation  of  the  light  irradiance  /, 
and  nonlinear  phase  shift  <p,  are: 


dl(z,t) 

dz 


al(z,t)-a2l2(z,t)  ^^-L^-  =  kn2I(z,t') 

and  dz 


(2) 


where  z  is  the  propagation  distance,  a  is  the  linear  absorption  coefficient,  ai  is  the  two-photon 
absorption  (2PA)  coefficient,  k  is  the  wave  number  and  m  is  the  nonlinear  refractive  index  from 
n=no+n2 1.  Here,  in  MKS  units,  which  are  used  throughout  this  chapter  unless  noted  otherwise, 


— Imj“  «,(«;«))  = -^4  Re 
"»  and  2"»  ,  (3) 

where  ^<3>  =  -%!;3)  (co;-a>,a>,a>)  f°r  a  single  linear  polarized  input  frequency  of  co,  where 

eff  2  '  ’xxxx 

E{z,t)=Ea  cos {kz-cot)  ,  and  Em(z)  is  the  slowly  varying  field  amplitude  with  irradiance, 

/  =  ^cs0n0 \Ea |2  .[9]  We  must  mention  here  that  this  definition  of  <22  in  terms  of  Im/13*  assumes 

that  we  are  not  within  a  region  of  linear  loss  where  other  effects,  including  absorption  saturation 
would  contribute.  Equations  3  assume  that  the  linear  absorption  <z«2cono/c,  otherwise  a2 
becomes  dependent  on  both  the  real  and  imaginary  parts  of  x<3),  as  does  m.  This  condition  is 
normally  satisfied  except  very  near  linear  resonance.  [10]  For  example,  in  the  visible  for  an  index 
of  ~1.5  this  becomes  a«  2xl05  cm"1;  however,  for  metallic  nanoparticles,  etc.,  this  can  be 
important. 

The  solution  of  Eq.  2  for  the  irradiance  propagation  is 


e-a,(\-R)2I(()j) 

\  +  a2(\-R)I(0j)L 


with  Lt^{\-QaL)/a  where  L  is  the  sample  length,  1(0, t)  is  the  irradiance  incident  on  the  sample, 
and  R  is  the  surface  reflectance.  If  the  thin  sample  approximation  is  not  valid,  amplitude  and  phase 
become  intertwined  with  propagation  within  the  sample  and  the  full  wave  equation  must  be 
solved.[ll]  In  the  following  sections  we  discuss  specific  methods  for  measuring  the  2PA 
coefficient  along  with  some  complementary  methods  for  determining  m.  NLR  always  occurs 
when  ce  is  being  measured,  and  sometimes  it  can  be  important.  In  addition,  because  2PA  and 
bound-electronic  m  are  related  through  Kramers-Kronig  relations  via  causality,  careful  and 
complete  measurements  of  nonlinear  spectra  can  be  used  to  derive  spectra/dispersion  of  the  other. 
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One  should  note  that  these  Kramers-Kronig  relations  are  strictly  only  valid  for  the  frequency 
nondegenerate  forms  of  the  nonlinearities.  [9] 

When  making  measurements  of  organic  dyes,  solutions  are  often  used  since  growing  crystals  can 
be  very  research  intensive  and  making  good  optical  quality  solid-state  films  can  also  be  difficult. 
In  such  measurements  nonlinear  cross  sections  are  the  more  meaningful  parameters  to  use.  These 
cross  sections  are  usually  given  in  GM  [12]  units  (1  GM  =  lxlO50  cm4s).  For  nonlinear  absorption 
the  cross  section  in  terms  of  the  2PA  coefficient  is 


a2hco  k0n2hco 

,  and  r=  N 


(5) 


where  we  have  defined  the  NLR  coefficient  in  an  analogous  fashion,  with  N  the  molecular  density. 
The  nonlinear  coefficient  per  molecule  may  also  be  defined  in  terms  of  the  second 
hyperpolarizability  y  which  is  defined  by: 


r 


(3) 


Nf4 


(6) 


where  /  =  (n2  +  2)/3  is  the  local  field  correction  factor.[13]  The  molecular  second 


Z(3) 

hyperpolarizability  is  usually  given  in  esu  as  Yesu  -  — ‘  4  .  For  noninteracting  solute  and  solvent 

CgsJ 

molecules,  assuming  a  nonabsorbing  solvent,  the  nonlinear  susceptibility  of  the  solution  is  given 
by: 


Z 


(3) 


£n 


(N,^r,.lm+NRtyf  +{Nlmrf 


1/2 


(V) 


where  the  index  used  in  the  local  field  correction  factor  is  the  index  of  the  solution,  i.e.  solvent 
plus  solute,  and  the  y's  are  the  orientationally  averaged  molecular  hyperpolarizabilities.  If  Woivcnt 
Ysoivent  =  -Absolute  Ysoiute  (the  contribution  of  the  real  parts  due  to  solute  and  solvent  cancel  each  other), 
the  absolute  value  of  xl3  )  of  the  solution  has  a  minimum  value  given  by:  [14] 


y(3) 

A/  solution 


V4NsoluteImysolute  =  \lmz: 


•(3) 

solution 


(8) 


This  allows  a  separation  and  determination  of  the  Re  and  Im  parts  of  the  nonlinear  susceptibilities 
by  varying  the  concentration  to  find  the  minimum. 
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Figure  1.  Left,  Schematic  energy-level  diagram  showing  two-photon  absorption  from 
ground  to  excited  state  with  equal  photons  (solid  green  lines)  of  energy  ftco,  or  unequal 
photons  (dotted  blue  and  red  lines)  of  energies  hm  and  tmi  where  2h($=tm\  +  ft  002.  The 
dashed  horizontal  lines  indicate  what  are  commonly  referred  to  as  virtual  states.  Middle, 
resonant  2-step  excitation  of  the  upper  state,  u,  with  linear  absorption  cross  sections 
between  the  states  of  agc  and  aeu  respectively.  The  rectangles  are  meant  to  indicate 
rovibronic  transition  absorption  bands.  Right,  two-photon  absorption  to  a  band,  e, 
followed  by  excited-state  absorption  to  band  u. 

2.2  Excited-State  Absorption  and  Refraction 

The  above  description  of  2PA  and  m  is  strictly  for  bound-electronic  nonlinearities  which  are 
ultrafast  (sub-femtosecond).  There  are  many  other  nonlinear  responses  which  are  pulsewidth 
dependent,  and  we  describe  just  a  few  of  these  in  the  following.  Many  of  these,  such  as  excited- 
state  nonlinearities,  involve  resonant,  i.e.  real  absorption  processes. [15]  Organic  dyes  often 
display  excited-state  nonlinearities  both  in  absorption,  ESA,  and  in  refraction,  ESR.  These  effects 
are  due  to  the  redistribution  of  level/state  populations.  Figure  1  shows  ‘instantaneous’  2PA  (left) 
where  the  intermediate  off  resonance  state  is  labelled  ‘e’  and  the  final  upper  state  is  ‘u’.  Transitions 
are  1 -photon  allowed  from  g  to  e  and  from  e  to  u  but  not  allowed  by  parity  from  g  to  u.  It  is 
instructive  to  look  at  the  specific  case  of  one-photon  induced  ESA  in  the  case  where  saturation  for 
the  g  to  e  transition  is  relatively  small. [15]  Figure  1  (middle)  shows  resonant  2PA,  i.e.,  the 
intermediate  state  resonantly  absorbs  /too  and  is  populated.  The  equations  describing  this 
phenomenon  in  terms  of  absorption  cross  sections,  a;  are: 


dl  Ar  r  Ar  _  dNe  cr  N I 


(9) 
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where  Ng  is  the  ground-state  population  density  and  Nc  the  population  density  of  excited-state  e  in 
Fig.  1  middle,  and  the  cross  sections  are  between  the  labelled  1 -photon  allowed  transitions.  As 
the  absorption  has  not  been  written  in  terms  of  a  population  difference  we  have  included  saturation 
of  the  linear  absorption  only  as  a  ground-state  depletion  to  low  order,  and  we  have  assumed  that 
once  the  excited  state  has  been  populated,  it  does  not  decay  back  to  the  ground  state  within  the 
duration  of  the  pulse.  This  is  often  a  good  approximation  since  the  initial  absorption  in  molecules 
is  usually  into  a  rovibronic  band  which  rapidly  relaxes  to  a  lower  energy  from  which  ESA  occurs, 
and  this  approximation  is  especially  good  when  considering  triplet  states  (not  shown). 


As  2PA  and  ESA  are  typically  measured  using  short  laser  pulses,  for  which  the  measured  quantities 
are  the  transmitted  pulse  energy  or  fluence  (energy  per  unit  area),  we  need  to  integrate  Eqs.  9  over 
time.  Under  our  approximations  these  two  equations  can  be  integrated  over  time  to  give  a  new 
propagation  equation: 


dF_ 

dz 


=  — <t  N  F  -  <7  N 

ge  g  ge  g 


ItlCO 


(10) 


where  F  is  the  time  integrated  irradiance  or  fluence,  F  =  7(t)dt.[15,  16]  Equation  10  is  good 

to  second  order  in  F.  The  last  term,  proportional  to  the  square  of  crge,  comes  from  ground-state 
depletion  which  can  sometimes  be  ignored  as  is  true  for  the  experimental  case  shown  in  Sec.  3.1.3. 
This  equation  is  analogous  to  the  equation  for  2PA  except  that  the  fluence  has  replaced  the 
irradiance  and  ai  is  replaced  by  products  of  one-photon  cross  sections.  When  performing 
perturbation  theory  to  calculate  two-photon  absorption,  it  is  the  product  of  one-photon  absorption 
matrix  elements  (proportional  to  the  linear  absorption  cross  sections)  that  enters  the  calculation. 
Thus  this  equation  shows  the  transition  from  nonresonant  2PA  to  resonant,  two-step,  nonlinear 
absorption.  There  is  a  limited  applicability  of  this  solution  as,  for  example,  when  the 
irradiance/ fluence  is  increased,  ground-state  depletion  becomes  important.  [16, 17]  Looking  at  the 
similarity  between  Eq.  10  and  Eq.  2,  it  is  not  surprising  that  this  ESA  has  been  mistaken  for  2PA 
in  some  experiments.  We  show  an  example  in  Sec.  3.1.3.  It  is  also  the  case  that  some  authors 
define  an  effective  2PA  coefficient  to  describe  this  ESA  process;  however,  this  coefficient  is 
pulsewidth  dependent  while  the  true  2PA  coefficient  is  independent  of  pulsewidth. 

Associated  with  this  ESA  is  ESR  since  a  change  in  absorption  also  results  in  a  change  in  index, 
both  due  to  the  redistribution  of  populations.  The  sign  of  the  NLR  will  depend  on  the  relative 
spectral  positions  of  the  resonances,  i.e.  above  resonance  or  below.  In  this  case  the  nonlinear 
phase  shift  is  given  by 

=kn2I(z,t)  +  (creuR  +  cr  )Ne(z,t)  ,  (11) 

az 

where  crgeR  and  <7eUR  are  the  nonlinear  refractive  cross  sections  for  the  ground  and  excited  state 
respectively  and  we  have  kept  the  instantaneous  NLR.  Here  the  integration  over  time  involves  an 
averaging  of  the  overall  phase  shift  which  results  in  a  similar  integral  as  for  ESA  which  does  not 
depend  on  pulse  shape,  giving  an  averaged  phase  shift  of 'A  the  phase  shift  at  the  peak  of  the  pulse, 
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dz 


{^geR  +  ^exR 


00 

f dtl(z,t)\‘  dt'I(z,t') 

\  a  j  J-co _ =  U 

/J-  pco  <-\V 


ESR 


1*00 

dtl(zj) 

J-00 


2~?geR+°e*R>hw 


(12) 


On  the  other  hand,  the  instantaneous  term  gives  an  average  phase  shift  that  depends  on  pulse  shape 
and  is  1/V2  for  Gaussian  temporal  pulses,  i.e., 


d{(/){z,t))t 


dz 


■  kn 


.L  (z?  0  kn2 


n2 


f  dtl{z^t) 

J— oo 


F 


(13) 


Measurements  of  2PA  in  semiconductors  and  many  organic  dyes  are  also  accompanied  by  ESA 
but  where  the  excited  state  is  populated  by  2PA.  In  this  case,  the  irradiance  propagation  equation 
becomes, 


=  _ al (z,  t)  -  a2I2  (z,  t)  -  <JexNeI ,  with 
dz 


dN e_ 
dt 


a21 2 
2  Tuo 


(14) 


as  discussed  in  more  detail  in  the  experimental  section  where  data  is  shown. 

This  is  a  higher-order  process  i.e.  2PA  followed  by  linear  absorption  (Iin^1*),  or  2PA 

followed  by  linear  refraction  (Re^1*)  as  indicated  in  Fig.  1,  right.  Thus  these  are  effectively  fifth- 
order  nonlinearities,  so  that  this  ESA  can  easily  be  confused  with  three-photon  absorption  (3PA), 
and  many  authors  have  defined  effective  3PA  coefficients,  but  again  these  will  be  pulsewidth 
dependent  by  necessity. [18]  These  nonlinearities  are  discussed  in  more  detail  in  the  sections 
describing  specific  experimental  arrangements,  e.g.,  Sec.  3.1.3  and  4.2. 

The  similarities  of  the  propagation  equations  for  ESA  and  2PA  as  well  as  similar  issues  for  higher- 
order  nonlinearities  present  challenges  for  characterization  techniques  to  uniquely  determine  the 
underlying  physical  processes.  Unique  determination  often  requires  multiple  methods  and 
parametric  studies,  e.g.  pulsewidth  dependence.  In  the  following  sections  we  discuss  several 
techniques  for  measuring  these  nonlinearities,  the  limitations  of  these  methods,  and  possible 
combinations  of  methodologies  to  eliminate  ambiguities. 

3.  METHODS  FOR  MEASUREMENTS  OF  NLA  AND  NLR 


Measurements  of  NLA  and  NLR  fall  into  two  broad  categories:  direct  methods,  where  we  measure 
the  self-induced  change  in  transmission  of  a  beam,  and  indirect  methods,  where  an  excitation  is 
induced  by  a  beam,  to  be  sensed  by  some  other  means  (transmittance  of  a  probe,  heating  of 
materials,  fluorescence,  etc.)[19] 

3.1  Direct  Methods: 


3.1.1  Nonlinear  Transmission: 

We  start  with  perhaps  the  most  straightforward  and  direct  measurement  of  transmission  as  a 
function  of  irradiance.  From  Eq.  4,  we  find 
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(15) 


i  _  mt) 

T  I(L,t ) 


-^[1  +  ^(0,  t){\-R)Leff] 
(l-K) 


where  T  is  the  transmittance. [20,  21]  The  equation  is  written  in  this  way  to  point  out  a  simple 
method  for  determining  the  2PA  coefficient  by  measuring  the  transmittance  and  plotting  inverse 
transmittance  versus  irradiance.  The  result  should  be  a  straight  line  whose  slope  determines  ai  and 
intercept  determines  a.  Of  course,  using  Gaussian  spatial  and  temporal  profiles  (or  other  shapes) 
requires  integration  over  space  and  time  These  integrals  reduce  the  slope  for  Gaussian  shapes  by 
a  factor  of  2V2  (root  2  for  each  dimension  of  space  and  time),  and  for  higher  irradiance  produces 
a  downward  curvature  due  to  spatial  and  temporal  beam  reshaping,  thus  deviating  from  a  simple 
straight  line. [22] 


Similarly  for  3PA 


1 

-im 

2_  f 

~M)T 

\  +  aj‘aui\-Rf 


''l  +  e-MY' 


a 


J) 


(16) 


where  it  is  assumed  that  there  is  no  2PA.  Figure  2  shows  a  plot  of  inverse  transmittance  for  a  two- 
photon  absorbing  organic  compound  in  thin- film  form  using  femtosecond  810  nm  pulses. [23] 

Only  the  low  irradiance  inputs  are  plotted  in  Fig.  2  since  at  higher  inputs  2-photon  excited  states 
linearly  absorb  leading  to  a  higher-order  nonlinear  response.  The  slight  upward  turn  as  opposed  to 
a  slight  downward  curvature  predicted  due  to  the  beam  broadening  by  nonlinear  absorption  is  due 
to  this  higher  order  process  of  2-photon  excited-state  absorption.  The  straight  line  in  Fig.  2  gives 
a  slope  which  yields  a  value  of  oa  =  30  cm/GW  which  agrees  with  the  values  reported  in  Ref.  [23], 

While  the  experimental  procedure  appears  to  be  simple  and  straightforward,  one  must  beware  of 
difficulties  that  result  from  the  presence  of  nonlinear  refraction  in  these  experiments.  If  this 
nonlinear  refraction  is  defocusing,  it  is  quite  easy  for  some  of  the  beam  to  miss  the  detector  at 
higher  irradiances  giving  the  appearance  of  a  larger  nonlinear  loss.  It  is  therefore  not  surprising 
that  most  inconsistencies  in  the  literature  involve  reported  values  for  the  nonlinear  coefficients 
being  too  large.  It  is  also  important  that  the  detector  response  is  spatially  uniform.  Additionally, 
as  for  any  nonlinear  measurement,  the  irradiance  needs  to  be  carefully  calibrated  necessitating 
measurements  of  the  spatial  profile  as  well  as  the  temporal  profile.  Gaussian  spatial  profile  beams 
are  desirable  since  their  propagation  is  particularly  simple. 
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Figure  2.  Data  on  a  25  micron  thick  film  of  MEH-PPV;  DOP  (50:50  melt-processed  blend 
where  DOP  is  di-octylphthalate,  a  plasticizer)  as  measured  in  Ref.  [23],  data  provided  by 
J.  Hales,  and  replotted  as  inverse  transmission  versus  irradiance.  The  upper  left  is  MEH- 
PPV,  and  is  lower  right  is  DOP.  The  wavelength  used  was  810nm  with  110  fs  (FWHM) 
pulses. 


The  literature  is  filled  with  reports  of  high  2PA  coefficients  that  are  the  result  of  either  missing  the 
beam  due  to  nonlinear  refraction  or  the  use  of  pulses  that  have  rapid  temporal  modulation  that  was 
not  noticed.  [22,  24]  Similar  problems  can  occur  if  the  spatial  profile  is  not  smoothly  varying  at 
the  sample.  These  comments  can  be  applied  to  nearly  all  of  the  measurement  techniques  described 
in  this  chapter;  however,  there  are  ways  to  mitigate  such  problems,  and  many  of  these  involve 
using  known  nonlinear  samples  as  references  or  test  samples. 

Another  comment  about  single  beam  direct  transmission  measurements  that  is  self-evident  is  that 
they  are  sensitive  to  all  nonlinear  absorption  processes  and  they  do  not  differentiate  between  the 
different  mechanisms  that  give  rise  to  them;  thus,  other  techniques  must  usually  be  used  to 
determine  the  process  or  processes  involved. 

3.1.2  Z-scan: 

Introduced  in  1989,  [25]  the  original  Z-scan  consists  of  translating  a  thin  sample  through  the  focus 
of  a  laser  beam  while  monitoring  the  light  transmitted  through  an  aperture.  [26]  This  gives  a 
measure  of  the  nonlinear  refraction  if  there  is  no  nonlinear  absorption  present.  Gaussian  spatial 
beams  are  preferred  since  modeling  the  results  is  greatly  simplified,  but  not  necessary.  For 
measuring  ultrafast  nonlinearities  pulsed  sources  are  required  and  the  signal  is  proportional  to  the 
temporally  integrated  nonlinearity.  Figure  3  (left)  shows  the  usual  Z-scan  setup  for  simultaneously 
performing  open  and  closed  aperture  Z-scans,  while  (right)  shows  the  configuration  for  a  closed- 
aperture  Z-scan  where  a  reference  arm  is  utilized  to  increase  the  signal-to-noise  ratio,  S/N.[27] 
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When  nonlinear  absorption  is  to  be  measured  care  is  taken  to  open  the  aperture  and  collect  all  of 
the  transmitted  light.  Typical  signal  profiles  for  Z-scans  are  shown  in  Fig.  4. 


Figure  3.  (Left)  Z-scan  setup  allowing  for  simultaneous  open  and  closed  aperture  Z- 
scans,  (Right)  Closed  aperture  Z-Scan  with  identical  reference  arm  increasing  the 
signal-to-noise  ratio. [27] 


Z /% 


Figure  4.  Calculated  Z-scan  signals  using  open,  closed  and  divided  signals  as  described 
for  third-order  nonlinearities  (2PA  and  self-defocusing)  with  the  parameters  shown  on 
the  figure  using  the  analysis  described  in  Sec.  3. 1.2.1  The  separation  of  the  transmission 
“peak”  and  “valley”  for  the  divided  signal  is  typically  labelled  A Zpv,  while  the  difference 
between  peak  and  valley  transmittances  is  commonly  labelled  A Tpv. 
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A  beam  splitter  is  often  placed  after  the  sample  shown  in  Fig.  3  to  split  the  transmitted  beam  onto 
two  detectors.  The  first  detector  measures  the  total  transmitted  light,  and  is  referred  to  as  an  open- 
aperture  Z-scan.  The  second  beam  passes  through  the  aperture  in  the  far-field  which  typically 
transmits  ~40%  of  the  light  in  the  linear  regime  which  is  collected  onto  a  second  detector.  This  is 
referred  to  as  a  closed-aperture  Z-scan,  i.e.  partially  closed.  The  open-aperture  Z-scan  is  sensitive 
only  to  nonlinear  absorption  as  long  as  the  sample  is  thin  as  previously  defined  and  the  curve  can 
be  fit  with  a  parameter  describing  the  NLA  (assuming  there  is  a  single  NLA  process!).  The  closed- 
aperture  Z-scan  is  affected  by  both  nonlinear  refraction  and  absorption  and  given  that  the  NLA  has 
previously  been  fit,  can  be  fit  with  a  single  parameter  for  the  NLR  (again  under  the  assumption 
that  there  is  a  single  NLR  process!).  However,  if  the  closed-aperture  signal  is  divided  by  the  open- 
aperture  signal  the  resulting  curve  is  nearly  identical  to  the  signal  that  would  have  been  obtained 
if  the  nonlinear  absorption  were  absent  -  see  Fig.  5.  This  can  greatly  simplify  data  analysis 
although  the  curve  may  also  be  fit  with  a  2-parameter  fitting  procedure  for  the  nonlinear  absorption 
and  nonlinear  refraction  coefficients  for  third-order  response  materials,  but  with  the  usual 
increased  errors  associated  with  a  2-parameter  fit.  The  popularity  of  the  Z-scan  technique  is  due 
in  large  part  to  the  simplicity  of  separating  absorptive  and  refractive  nonlinearities;  however,  it  is 
sensitive  to  ALL  NLA  and  NLR  mechanisms. 


Figure  5.  Z-scans,  Open  (black  squares),  Closed  (red  circles),  Divided  (blue  triangles)  for 
the  squaraine  molecule  shown  along  with  fits  using  only  ultrafast  m  and  «2.[28] 

Shown  in  Figure  5  is  an  example  of  Z-scans  (actually  Dual-Arm  Z-scans  described  in  Sec. 

3. 1.2.2)  using  femtosecond  pulses  where  only  2PA  and  m  are  observed.  The  nonlinear  refraction 
is  from  the  solute  only  where  the  NLR  from  the  solvent  has  been  subtracted  as  described  in  Sec. 
3.1 .2.2.  In  the  usual  Z-scan,  in  order  to  determine  the  m  of  the  solute  the  NLR  of  the  solvent 
must  be  subtracted  by  taking  two  sequential  Z-scans,  one  with  only  the  solvent  and  the  other 
with  the  solution. 
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3. 1.2.1  Simple  Z-scan  analysis: 

A  simple  empirical  relation  exists  between  the  peak  to  valley  change  in  transmission  measured  in 
a  Z-scan  with  Gaussian  beams  for  the  closed  aperture  Z-scan  signal  and  the  nonlinearly-induced 
phase  distortion:  [26] 

A7;,=  0.406(1  -S)“|A®„|,  (17) 


where  S  is  the  linear  transmittance  of  the  aperture  placed  in  the  far  field,  and, 

AO0(t)  =  ^n2I0(t)Leff, 


(18) 


where  7o(t)  is  the  an  axis  irradiance  at  focus  and  the  effective  length,  Lefi={\-Q'aL)l  a.  This  empirical 
relation  is  accurate  for  both  cases  without  NLA  and  with  NLA  after  division  by  the  open-aperture 
Z-scan  data.  Since  the  temporally  averaged  transmittance  is  measured  in  a  Z-scan,  the  pulse  shape 
needs  to  be  included.  This  is  relatively  straightforward  given  the  linear  relation  between 
transmittance  change  and  phase  distortion  of  Eq.  17.  If  the  irradiance  dependent  pulse  shape  in 
time  is  described  by  the  function  f( t),  then  the  temporal  integral  given  by 


At  =  f”  f2(t)dt  r  f(t)dt,  reduces  the  measured  A Tpv  by  Ax.  This  has  the  value  of  1/V2  for 

J— OO  /  J— 00 


Gaussian  pulses  for  nonlinearities  much  faster  than  the  pulsewidth  used.  For  pulses  much  shorter 
than  the  nonlinear  response,  e.g.  excited-state  nonlinearities,  Ax=l/2  independent  of  the  pulse 
shape.  In  addition,  the  separation  in  Z  between  the  position  of  the  peak  and  valley  is  given  by: 


|AZpv|*1.71Z0,  (19) 

where  Zo=»mi’o2/A,  with  wo  the  beam  HWl/e2M  in  irradiance.  This  allows  an  independent 
determination  of  the  focused  spot  size  if  the  sample  is  known  to  exhibit  a  third-order  nonlinear 
response. 

For  the  open  aperture  Z-scan,  the  peak  energy  transmittance  change  is  given  by: 


AT(Z)  *  - 


<h  1 
2V2  [1  +  Z2/Z2]’ 


(20) 


where  qo  =a2lo(0)Leff  (|go|«l).  Here  no  propagation  to  the  far  field  is  required  in  the  calculation 
since  all  energy  is  detected.  The  irradiance  at  the  exit  of  the  sample  7e  is: 


7e(Z,r,t) 


7(Z,r,t)exp~aL 
1  +  q(Z,r,t) 


(21) 


while  the  phase  change  at  the  exit  is: 


KYI 

A(j){Z,  r,  t )  =  — -  ln[l  +  q(Z,  r,  i)] , 
a2 


(22) 
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where  q(Z,r,t)=a2l(Z,r,i)Leff.  This  phase  along  with  the  irradiance  distribution  must  be  propagated 
to  the  detector  in  the  far  field  to  model  the  closed-aperture  Z-scan  data.  There  are  multiple  ways 
to  perform  this  propagation  and  two  of  these  are  described  in  Ref.  [26].  Thus,  while  the  empirical 
results  above  are  useful,  the  full  propagation  results  can  be  fit  for  all  Z  which  is  good,  for  example, 
for  checking  the  order  of  the  nonlinearity.  For  example,  a  3PA  signal  will  be  narrower  in  Z  than 
a  2PA  signal  and  |AZpv|  =1.5Zo  for  a  fifth-order  nonlinearity. 

Combining  Eqs.  21  and  22  we  obtain  the  complex  field  at  the  exit  surface  of  the  sample  to  be 

EXz,r,t)  =  E{Z,r,t)e^\\  +  q{Z,rj)f^la^n)  ,  (23) 


while  the  field  at  the  aperture  is  found  by  performing  a  Fresnel  diffraction  calculation: 


ATT  7  TTY  i*oo 

E.  (Z, r,  t,  d)  = - exp  -  [  r'dr'  E  (Z, r' ,  t  -  d'/c)  exp 

ad'  I  *.d'  Jo  e 


f  inr2 


A  d 


f i7rr'2)] 

^2;rrr'^ 

Ud’J  ° 

Ud'  J 

(24) 


where  d,=d-Z  is  the  distance  from  the  sample  to  the  aperture  plane.  The  measured  quantity  is  the 
pulse  energy  or  average  power  transmitted  through  the  far-field  aperture  having  a  radius  of  ra.  The 
normalized  transmittance  is  then  obtained  as: 


=  £>f|£.(Z,r,t,d)|2«fr 

1  ’  u 

where  U  is  the  same  as  the  numerator  but  in  the  linear  regime  (i.e.  for  A0  =  0  ).  In  the  case  of  an 
EZ-scan,  the  limits  of  the  spatial  integral  in  Eq.  25  must  be  replaced  by  rd  to  oo  where  rj  is  the 
radius  of  the  obscuration  disk.[29,  30]  It  is  generally  more  convenient  to  represent  the  aperture 
(or  disk)  size  by  the  normalized  transmittance  (or  rejection)  S  in  the  linear  regime. 


The  formalism  thus  far  presented  is  generally  applicable  to  any  radially  symmetric  beam.  Here, 
however,  we  assume  a  TEMo.o  Gaussian  distribution  for  the  incident  beam  as  given  by: 


wn 


E(Z,  r,  t)  =  E()  (t)  exp 
w(Z) 


■  +  l  - 


nr 


w/(Z)  AR(Z) 


+  i(f> 


(26) 


where  w(Z)=wo(l+Z2/Zo2) 1/2  and  R{Z)=  Z+Zo2/Z.  The  radially  invariant  phase  terms,  contained  in 
(/>,  are  immaterial  to  our  calculations  and  hence  will  be  ignored. 


The  integral  in  Eq.  24  can  be  analytically  evaluated  if  we  assume  that  |g|<l  (in  Eq.  23)  and  then 
perform  a  binomial  series  expansion  of  Ee  in  powers  of  q.  Recalling  that  q  oc  /  oc  exp(-r2/>v2),  this 
expansion  effectively  decomposes  Ee  into  a  sum  of  Gaussian  beams  with  varying  beam  parameters. 
This  method  of  beam  propagation  known  as  Gaussian  decomposition  was  first  given  by  Weaire  et. 
al.[31]  Following  the  expansion,  we  obtain: 

oo 

Ee  =  E{Z,r,t)e~aLI2YJFm  exp(2 mr2/w2(Z)),  (27) 

m= 0 
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where  Fm ,  the  factor  containing  the  nonlinear  optical  coefficients,  is  given  by: 


m\ 


j= i 


1  +  i 


n  aj3 


J  2  J  2m2 


(28) 


Remarkably  this  single  beam  technique  has  a  demonstrated  sensitivity  to  induced  phase  distortion 
of  ~/l/103  while  a  simple  variation,  the  “EZ-scan”,  which  replaces  the  aperture  with  a  stop,  has  a 
demonstrated  sensitivity  to  Z/104.[29]  There  have  been  multiple  other  modifications  of  the  basic 
Z-scan  reported  in  the  literature  [27,  32-36].  This  interferometric-like  sensitivity  comes  about  by 
the  sample  serving  as  a  phase  mask,  which  upon  propagation  to  the  far  field  is  transformed  into  a 
spatial  amplitude  redistribution,  i.e.  via  diffraction.  That  is,  diffraction  is  an  interference 
phenomenon  between  the  wings  of  the  beam  which  are  weak  and  propagate  nearly  linearly  and  the 
center  of  the  beam  which  undergoes  a  nonlinear  phase  shift,  either  advanced  as  for  self-defocussing 
(«2<0),  or  delayed  as  for  self-focusing  («2>0).  Thus,  the  Z-scan  serves  as  a  single  beam 
interference  technique  without  the  alignment  complexities  of  interferometers. 

3. 1.2.2  Dual  Arm  Z-scan 


Shown  in  Figure  5  is  an  example  of  Z-scans  where  only  2PA  and  m  are  observed.  Most,  solvents 
do  not  exhibit  2PA  in  the  visible  or  near  infrared,  so  any  observed  2PA  can  be  attributed  to  the 
solute.  However,  the  same  is  not  true  of  NLR,  and  any  measured  nonlinear  refraction  comes  from 
a  combination  of  the  NLR  of  the  solvent  and  solute.  In  order  to  determine  the  m  of  the  solute,  the 
NLR  of  the  solvent  must  be  subtracted.  Since  sometimes  the  solubility  of  molecules  is  limited,  the 
solvent  n2  can  overwhelm  the  solute  NLR  and  make  its  determination  problematic.  One  method  to 
alleviate  this  problem  is  the  "Dual- Arm”  (DA)  Z-scan  [28]  described  below. 

In  the  DA  Z-scan  (experiment  shown  in  Fig.  6),  a  second  Z-scan  arm  is  arranged  to  have  an 
identical  optical  path  to  that  of  the  first  arm.  When  initially  aligned  with  identical  cuvettes  and  the 
same  solvent  in  each  arm,  the  signal,  which  is  the  subtraction  of  the  signals  from  each  arm,  can  be 
nulled.  When  the  solute  is  added  to  one  arm,  the  resultant  difference  signal  is  then  only  due  to  the 
solute.  In  practice  we  find  that  this  increases  the  signal -to-noise  for  NLR  by  a  factor  of  ~10x  (see 
Figure  7  of  a  squaraine  molecule  showing  the  increased  S/N).  The  S/N  increase  for  NLA  is  less 
dramatic.  This  increase  in  S/N  for  NLR  comes  about  because  much  of  the  noise  in  each  arm  is 
correlated  on  each  laser  firing.  This  noise  is  effectively  subtracted  on  each  laser  shot.  For  example, 
since  closed-aperture  Z-scans  are  sensitive  to  beam  pointing  instabilities,  the  dual  arms  provide 
an  excellent  way  to  reduce  their  effect  since  both  arms  move  synchronously  thus  giving  the  same 
change  in  signal  from  this  beam  movement.  This  technique  also  works  well  for  thin  film  samples 
on  substrates  to  subtract  the  substrate  signal  and  correlated  noise.[37]  Figure  7  shows  the 
effectiveness  of  using  this  method  for  a  squaraine  dye  in  toluene  showing  both  NLA  and  NLR. 
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Figure  6.  Schematic  of  dual-arm  Z-scan.  The  items  labeled  CA  and  OA  represent  the 
closed  aperture  and  open  aperture  detectors  for  each  arm,  respectively.  The  reference 
beam  used  for  energy  monitoring  is  not  shown  (from  Ref.  [28]). 


Figure  7.  (a)  Sequential  CA  single-arm  Z-scans  of  the  solvent  toluene  (open  red  triangles) 
and  the  solution  of  SD-0  2405  in  toluene  (closed  black  squares)  at  695  nm  where  the 
concentration  C  =  47  pM,  the  pulse  energy  E  =  31  nJ  (Io  =  51  GW/cm2)  and  (b)  the 
subtraction  of  the  solvent  CA  signal  from  the  solution  CA  signal  (open  green  squares); 
Note  this  is  for  the  same  molecule  as  shown  in  Fig.  5.  (c)  Simultaneous  CA  dual-arm  Z- 
scans  of  the  solvent  toluene  (open  red  triangles)  and  the  solution  SD-0  2405  in  toluene 
(closed  black  squares)  at  695  nm  using  the  same  pulse  and  (d)  the  subtraction  of  the 
solvent  CA  signal  from  the  solution  CA  signal  after  low-energy  background  signal 
(LEB(Z))  subtraction  (open  green  squares)  and  corresponding  fit  of  both  2PA  and  NLR 
(solid  blue  line)  with  A</>o  =  -0. 16,  qo  =  0.077,  using  S  =  0.33. [28] 
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3.1.3  Determining  nonlinear  response  from  pulsewidth  dependence  of  Z-scans 

Z-scans  of  materials  exhibiting  one-photon  absorption  (1PA)  induced  ESA  (see  Eqs.  9-12)  cannot 
easily  be  distinguished  from  Z-scans  performed  on  materials  with  only  ultrafast  nonlinearities. 
One  way  to  distinguish  is  to  perform  the  Z-scan  at  two  different  pulsewidths.  If  the  Z-scans  are 
identical  for  the  same  irradiance  but  different  pulsewidths,  the  response  time  of  the  nonlinearity  is 
shorter  than  the  temporal  width  of  the  pulses  used.  For  very  short  pulses,  (e.g.  femtosecond)  this 
would  indicate  that  the  mechanism  is  2PA.  If  the  open  aperture  Z-scans  are  identical  for  the  same 
fluence,  then  they  may  be  due  to  ESA.  Similarly,  if  the  closed  aperture  Z-scans  are  the  same  for 
the  same  fluence  rather  than  irradiance,  then  the  nonlinear  refraction  may  be  due  to  the 
redistribution  of  the  populations. 

In  order  to  illustrate  this,  we  show  Z-scan  data  for  a  sample  of  chloro-aluminum  phthalocyanine 
(CAP)  in  Fig.  8  using  a  fixed  energy,  but  pulswidths  differing  by  a  factor  of  two.  [15]  Thus  the 
irradiance  changes  by  a  factor  of  two  while  the  signals  for  both  nonlinear  absorption  and  nonlinear 
refraction  are  nearly  identical.  Both  nonlinear  absorption  and  refraction  are  due  to  excited  states, 
rather  than  to  2PA  and  bound-electronic  n2. 


Figure  8.  Open  aperture  Z-scan  data  at  532  nm  using  (open  triangles)  30  ps  pulses 
(FWHM),  and  (open  squares)  62  ps  pulses  (FWHM)  on  Chloro-Aluminum 
Phthalocyanine  -  left,  and  closed  aperture  Z-scan  data  (divided  by  open  aperture 
data)  for  the  same  pulsewidths  -  right.  The  molecule  is  shown  at  the  right.  From 
Ref.  [15]. 

Creating  new  linear  absorption  resonances  by  creating  excited  states  changes  the  index,  and  the 
sign  depends  upon  which  side  of  the  resonance  the  light  is  -  and  additionally  removing  absorbers 
from  the  ground  state  also  changes  the  index  -  again  the  sign  depends  on  the  resonance  frequency 
and  these  contributions  add  to  determine  the  net  NLR.  These  can  be  determined  by  the  usual 
Kramers-Kronig  relations  for  the  materials  with  the  changed  absorption  profile,  i.e.  including  the 
new  ESA  or  saturated  absorption  profile.  It  is  important  to  note  that  this  is  not  an  m  since  it  will 
depend  on  the  density  of  excited  states  created,  which  is  energy  dependent. 
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When  the  excited  states  are  created  by  2PA,  as  described  in  Eq.  14,  the  NLA  can  mimic  fifth-order 
3PA  while  the  nonlinear  refraction  also  appears  as  a  fifth-order  NLR  signal.[18,  38]  Figure  9 
shows  a  molecule  where  ESA  and  excited-state  refraction  occur  due  to  2PA. 


Figure  9.  Open  (open  circles),  closed  (closed  circles)  and  divided  (closed  squares) 

Z-scan  data  of  the  organic  molecule  shown  on  the  right.  The  solid  lines  are  fits  to 
the  data.  From  Ref.  [38]. 

Here  the  equations  are  modified  in  a  straightforward  way  following  Eq.  14  but  ignoring  linear 
absorption.  The  irradiance  varies  as 

^p-  =  -a1l1{t)-caN,(  i)I(t)  (29) 

oz 


while  the  excited-state  generation  varies  as 

dNe  _a2I2(t) 
dt  2  hco 


(30) 


While  these  don't  have  a  simple  solution  as  in  the  1PA  generated  excited-state  case,  they  can  be 
numerically  evaluated  for  Z-scans  as  shown  in  the  fits  of  Fig.  9.  Again,  for  a  fixed  irradiance 
longer  pulses  contain  more  energy,  and  the  role  of  the  ESA  and  ESR  is  increased.  Some 
researchers  have  quoted  effective  3PA  coefficients  for  the  NLA  portion  of  these  equations  by 
temporally  integrating  the  excited-state  density.  For  example,  assuming  Gaussian  temporal  pulses 
of  HWl/eM  to  yields; 

=  f  l\d)df=^^\te~2{t'lto)2df  (31) 
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which  can  be  used  in  Eq.  29.  If  the  pulses  are  long  enough,  the  second  term  in  Eq.  29  dominates, 


I K  (X-,llt, 


and  since  Ne(t)  oc  I2,  the  open  aperture  Z-scan  signal  looks  like  3PA  with  Ne( co)  =  . 


21  0  l0 


2  2fuo 


which  when  used  in  Eq.  29  yields  an  73  dependence  for  the  spatial  derivative  of  the  irradiance.  This 
provides  a  quick  estimate  of  the  size  of  the  nonlinear  absorption.  However  it  also  points  out  how 
easily  this  cascaded  2PA+ESA  (i.e.  effect  can  be  misinterpreted  in  an  experiment  as 

3PA.[18]  The  nonlinear  refraction  from  creating  these  excited  states  also  looks  like  a  fifth-order 
"m  using  similar  math  (Im^^Re^1*);  however,  again  this  doesn't  give  the  essential  physics 
behind  the  NLR. 


3.1.4  White-Light-Continuum  Z-scan  (WLC  Z-scan): 

Performing  nonlinear  spectroscopy  using  Z-scan  can  be  time  consuming  since  most  broadly 
tunable  sources,  e.g.  optical  parametric  sources,  require  lengthy  characterization  each  time  they 
are  tuned  since  the  temporal,  and  in  particular,  the  spatial  profiles  are  not  reproducible.  Usually 
spatial  filtering  is  required  if  Gaussian  profile  pulses  are  to  be  used.  Again  this  is  not  a  necessary 
condition  for  performing  Z-scans  but  it  greatly  simplifies  analysis. [39-41] 

Nonlinear  spectroscopy  using  Z-scan  is  possible  as  shown  in  Figure  10,  where  the  2PA  spectrum 
of  a  squaraine  molecule  shown  in  the  figure  was  measured  by  performing  open-aperture  Z-scans 
at  many  wavelengths  using  an  optical  parametric  amplifier  (OP A).  Such  a  spectrum  can  take 
several  days  to  perform,  since  the  Z-scan  needs  to  be  realigned  for  each  wavelength,  as  the  OPA 
beam  alignment  changes  upon  tuning.  [42] 
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Figure  10.  Linear  absorption  (right  arbitrary  axis)  along  with  2PA  cross  section  (left 
axis)  determined  by  open  aperture  Z-scans  versus  photon  energy  for  the  molecule 
shown.  The  energy  corresponds  to  the  energy  of  final  excitation  for  both  linear  and  two- 
photon  absorption,  e.g.  the  2PA  peak  is  near  an  input  photon  energy  of  2.35/2  eV.[42] 

A  solution  that  allows  rapid  Z-scan  characterization  of  materials  over  a  broad  spectral  range  is  the 
white-light-continuum,  WLC,  Z-scan.  [3 9-41,  43]  Here  the  Z-scan  methodology  is  the  same  but 
we  replace  the  usual  'single'  frequency  source  in  the  Z-scan  with  a  spectrally  filtered  WLC. [44] 

The  apparatus  is  shown  in  Fig.  11.  The  WLC  is  produced  by  weak  focusing  in  a  long  ~1.5  M  cell 
filled  with  Kr  gas  at  a  pressure  of  -2  atm.  Using  Gaussian  spatial  profile,  femtosecond  excitation 
pulses  of  -lmJ  at  ~800nm  produces  a  clean  spatial  profile  pulse  of  spectral  content  from  <400nm 
to  >800nm,  i.e.  greater  than  1  octave. 

Kr  (2.4atm) 


Figure  11.  WLC  Z-scan  experimental  setup:  L:  lens;  M:  mirror;  WP:  half- wave-plate;  P: 
polarizer;  FW:  filter  wheel;  BS:  beam  splitter;  D:  detector;  A:  aperture;  S:  sample;  the 
dotted  M's  are  removable  mirrors  for  beam  characterization.  From  Ref.  [40] 
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The  entire  WLC  cannot  be  used  in  a  single  Z-scan  since  besides  the  degenerate  2PA  there  will  be 
strong  nondegenerate  2PA  and  these  two  processes  cannot  be  simply  separated.  The  simplest 
method  to  apply  is  to  spectrally  filter  the  WLC  prior  to  the  sample  and  perform  a  normal  single 
frequency  Z-scan.  Then  the  spectral  filter  can  be  changed  (e.g.  using  a  computer  controlled  rotating 
platform  holding  many  spectral  filters  as  shown  in  Fig.  11,  or  using  a  linear  variable  filter)  and 
another  Z-scan  performed. 

The  filters  must  be  spectrally  broad  enough  to  support  the  pulses  used,  in  our  case  ~  100  fs.  The 
difference  between  using  the  WLC  and,  for  example,  tuning  an  optical  parametric 
generator/amplifier,  OPG/A  is  in  the  length  of  time  it  takes  to  measure  a  nonlinear  spectrum  (or 
nonlinear  dispersion  curve).  As  it  turns  out,  once  characterized  in  terms  of  spatial  properties  at  all 
wavelengths  and  temporal  profiles  at  all  wavelengths,  the  WLC  is  stable  from  day  to  day  and  week 
to  week.  This  cannot  be  said  of  optical  parametric  sources.  We  have  found  that  OPG/A's  when 
tuned  change  their  spatial  profile  and  usually  spatial  filtering  is  needed  to  obtain  Gaussian  profile 
beams.  Figure  12  shows  the  characterization  of  a  WLC  demonstrating  that  near  Gaussian  profile 
beams  are  obtaining  at  all  the  wavelengths  over  an  octave  span.  New  research  is  extending  the 
useful  range  to  >  two  octaves  using  different  excitation  wavelengths.  Thus,  once  characterized, 
the  WLC  can  be  used  to  perform  nonlinear  spectroscopy  while  the  OPG/A  needs  to  be  re¬ 
characterized  after  each  tuning  cycle. 
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Figure  12.  Characterization  of  the  WLC.  Spatial  profiles  -  upper  left,  temporal  pulsewidth 
-  upper  right,  available  energy  after  narrow  band  filters  -  lower  left,  example  of  an 
autocorrelation  after  a  620nm  narrow  band  filter  (~10  nm  bandwidth)  -  lower  right. [40] 

This  WLC  Z-scan  has  the  potential  of  serving  as  the  nonlinear  equivalent  of  a  linear 
spectrophotometer.  When  automated,  a  researcher  can  insert  a  sample,  initiate  a  measurement 
and  then  after  some  time  (perhaps  <  1  hour  for  kHz  sources)  retrieve  the  degenerate  nonlinear 
absorption  spectrum;  however,  with  this  nonlinear  spectrophotometer  there  is  the  added  benefit 
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not  given  by  a  linear  spectrometer.  The  dispersion  of  the  nonlinear  refraction  is  measured 
simultaneously  with  the  closed-aperture  Z-scan  as  shown  in  Fig.  1 1 .  Using  this  configuration  the 
data  of  Fig.  13  were  taken. 


FkCO 


Figure  13.  2PA  cross  section,  right  axis,  and  nonlinear  refractive  cross  section,  left  axis 
for  the  molecule  shown  to  the  right.  Data  taken  using  the  WLC  Z-scan,  except  for 
wavelengths  longer  then  800nm  where  2-photon  fluorescence  is  used  (open  circles). 
From  Ref.  [40] 

3.1.5  Other  Variants  of  the  Z-scan  method: 

Determining  the  physical  mechanism  or  mechanisms  responsible  for  the  signal  is  seldom  revealed 
by  a  single  Z-scan  experiment  although  using  ultrashort  femtosecond  pulses  can  sometimes  limit 
the  significant  responses  to  bound-electronic  nonlinearities.  [45,  46]  See  Fig.  8  where  either  one 
of  the  Z-scans  could  have  been  used  to  extract  erroneous  values  for  a  2PA  coefficient  and  an  m. 
There  have  been  many  variations  of  the  Z-scan  introduced  over  the  past  years.  Many  of  these  have 
shown  the  potential  to  further  increase  the  sensitivity,  e.g.  EZ-scan  [29],  while  others  have 
introduced  multiple  beams  which  have  the  potential  to  allow  for  measuring  nondegenerate 
nonlinearities  [27,  35,  47].  However,  all  these  “improvements”  come  at  the  cost  of  complicating 
the  methodology,  and  the  simplicity  of  this  technique  is  the  likely  reason  it  became  so  popular. 

The  4f  coherent  imaging  method  [48]  and  many  other  techniques  are  closely  related  to  Z-scan. 
Here,  for  example,  a  beam  profile  is  focused  through  a  nonlinear  sample  and  imaged  in  the  linear 
regime.  Then  the  irradiance  is  increased  and  the  change  in  the  imaged  irradiance  distribution  is 
observed.  The  nonlinearly  induced  phase  distortion  can  be  calculated  and  the  nonlinear  index 
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determined.  Nonlinear  absorption  can  also  be  analyzed.  In  the  particular  reference  given  [48],  a 
circular  flat-top  beam  is  imaged  and  the  light  going  into  the  originally  dark  area  in  the  image  plane 
is  monitored.  When  fully  aligned  the  nonlinear  index  (and/or  absorption)  can  be  determined  in  a 
single  shot. 

As  previously  stated,  the  Z-scan  technique  is  sensitive  to  ALL  NLA  processes  and  ALL  NLR 
mechanisms.  Unraveling  the  various  physical  processes  that  contribute  to  these  nonlinearities 
cannot  be  unambiguously  determined  with  single  Z-scans. 

3.2  Indirect  Methods 

By  indirect  we  mean  that  a  material  is  excited  optically  to  be  sensed  by  some  other  means  than  the 
transmission  of  the  excitation  beam  itself.  These  other  means  could  include  measuring  the 
transmittance  of  a  probe,  heating  of  materials,  fluorescence,  etc.  For  example,  an  excitation  beam 
induces  excited  states  via  2PA  whose  presence  are  sensed  by  their  fluorescence  whose  intensity  is 
proportional  to  how  many  excited  states  were  produced  via  2PA.  We  give  examples  of  these 
indirect  methods  in  the  following. 

3.2.1  Excitation-Probe  Methods: 

Figure  14  shows  the  standard  excite-probe  experimental  setup  where  an  excitation  pulse  is  incident 
on  the  system  at  time  Zd  =  0  and  the  probe  pulse  examines  the  change  in  transmission  properties 
of  the  sample  at  later  times  Zd  >  0.  For  example,  see  articles  in  Ref.  [49].  The  excitation  and  probe 
pulses  could  be  at  the  same  (degenerate)  or  different  (nondegenerate)  wavelengths  but  in  order  that 
they  are  synchronized  in  time  they  must  be  derived  from  the  same  pulsed  laser  source.  The  time 
delay,  zd,  of  the  probe  is  controlled  by  a  delay  line  as  shown  in  Fig  14. 


Figure  14.  Excite-probe  experimental  setup  showing  the  optical  delay  line. 

Two-photon  absorption  can  easily  be  differentiated  from  ESA  using  excite-probe  techniques  since 
for  times  after  short  pulse  excitation  (pulsewidth  shorter  than  the  decay  time  of  the  system  r),  the 
ESA  will  be  present  after  the  excitation  pulse  is  over,  i.e.  for  r  greater  than  the  pulsewidth.  For 
example,  for  1PA  generated  excited  states  the  dynamics  are  described  by  a  modification  of  Eq.  9 
that  includes  simple  decay 

dN  a  N I  N 

_  ge  g  e  (32) 

dt  ha>  z 
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where  definitions  are  as  given  in  Sec.  2.1 ,  Ie  is  the  excitation  irradiance  and  x  is  the  decay  time.  It 
is  assumed  that  the  probe  irradiance  is  weak  enough  that  it  generates  no  significant  excited-state 
population.  The  transmission  change  is  sensed  as  a  function  of  time  delay  by  the  probe  pulse, 
according  to 

“jf  =  ~{°geNg  +  ^euNe)lp-  (33) 


Provided  the  excitation  and  probe  pulses  are  much  shorter  than  the  decay  time,  r,  integration  over 
time  to  find  the  probe  fluence  is  straightforward.  Defining  the  probe  fluence  transmittance  as, 
Tprobe(Jd)  =  0)  >  and  the  normalized  transmittance  as  T Ni(Jd)  =  ^ probed d) /T l, 

where  TL  is  the  linear  transmittance  of  the  probe,  we  find, 

TNL(j)  =  [fWL(0)](exp(-Td/T))  (34) 


where  TNL  (0)  = 


1  4 


{?eu  age)Ft 


2  ha) 


:(i 


—  p  CCqL 


>r 


is  the  probe  fluence  transmission  for  the  probe 


delay  immediately  after  the  excitation  pulse.  Here  we  have  assumed  that  a0  =  ogeNg  is  constant, 
i.e.,  small  depletion  of  the  ground  state.  Provided  the  changes  in  normalized  transmittance  are 
small,  i.e.  AT  =  1  —  TNL  «  1  then  Eq.  34  becomes 


ATnl(t)  =  ATNL(0)e-T«/\ 


(35) 


so  that  the  transmittance  change  can  be  simply  fit  with  an  exponential  decay  to  determine  the 
excited-state  lifetime.  Meantime  the  excited-state  absorption  cross  section,  aeu,  can  be  determined 
from  T jvl(0)-  Even  if  we  do  not  satisfy  the  conditions  for  Eq.  35,  we  may  still  fit  the  decay  with 
Eq.  34;  however,  provided  we  remain  in  this  low-excitation  regime  where  there  is  small  depletion 
of  the  ground  state  population,  we  do  not  need  to  know  the  molecular  density  to  obtain  oeu.  Figure 
15  shows  an  example  of  an  excite-probe  experiment  on  a  squaraine  molecule  (SD  2577)  at  a 
wavelength  of  532  nm  with  30  ps  pulses. [50] 


Delay  (ns) 


Figure  15.  Excite-probe  on  the  squaraine  molecule  shown  in  the  inset  (SD  2577)  at  a  wavelength 
of  532  nm,  showing  ESA  with  a  decay  time  of  0.9  ns.  Details  are  given  in  Ref.  [51]. 
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This  single  excite-probe  experiment  can  often  determine  the  singlet  manifold  molecular 
parameters,  but  it  is  usually  incapable  of  determining  the  triplet  yield,  cross  section,  and  lifetime 
which  are  often  important  for  organic  dyes.  This  can  be  done  using  a  double  excite-probe 
geometry  as  discussed  in  detail  in  Ref.  [52],  In  these  experiments  a  first  excitation  creates  a 
population  of  the  triplet  state  which  is  excited  again  after  allowing  sufficient  time  to  create  a  triplet 
population  and  then  probed  to  determine  these  molecular  constants. 

While  degenerate  excite-probe  methods  are  probably  best  for  fully  understanding  the  mechanisms 
involved  in  a  single -beam  experiment  such  as  Z-scan,  it  is  often  useful  to  perform  nondegenerate 
excite-probe  measurements,  where  the  excitation  and  probe  are  at  different  wavelengths. 
Experimentally,  this  has  the  advantage  that  the  weak  probe  can  be  easily  distinguished  from  the 
strong  excitation  beam  using  spectral  filters,  but  using  different  excitation  and  probe  wavelengths 
allows  other  measurements  to  be  performed.  This  allows,  for  example,  the  entire  excited-state 
absorption  spectrum  to  be  measured  by  tuning  the  probe  while  keeping  the  excitation  wavelength 
fixed. [53]  Also  the  excitation  beam  can  be  used  to  generate  excited  states  via  two-photon 
absorption  which  has  the  advantage  of  creating  a  more  uniform  excitation  through  the  depth  of  the 
sample.  Finally,  there  may  be  nondegenerate  2PA  involving  the  absorption  of  one  photon  from 
each  of  the  excitation  and  probe  beams.  In  this  case,  the  propagation  equations  for  the  probe  and 
excitation  beams  are, 


dip 

dz 


2tt2 (Wp,  U)e)/g/p  &ex(^p^Ng!p 


dig 

dz 


-a2(coeicoe)I% 


(36) 


where  a2((op;  coe)  is  the  nondegenerate  2PA  coefficient,  and  it  is  assumed  that  there  is  no  ESA 
experienced  by  the  excitation  beam.  In  the  event  that  there  is  degenerate  2PA  at  the  excitation 
wavelength,  the  excited-state  population  is  given  by, 


dNe  _  «2 ((QeWeVe 
dt  2hu> 


(37) 


As  illustrated  in  Fig.  1  (right),  the  effects  of  these  excited  states  may  be  seen  by  the  presence  of 
ESA  after  the  excitation  pulse,  and  since  the  excited-state  density  is  proportional  to  7e,  the  ESA 
experienced  by  the  probe  will  appear  as  a  higher-order  nonlinearity  than  the  nondegenerate  2PA, 
which  is  proportional  to  Ie.  With  femtosecond  pulses,  the  population  densities  generated  by  2PA 
are  usually  too  small  for  significant  ESA  to  be  observed.  However,  Fig.  16  shows  a  lead 
bis(ethynyl)porphyrin  polymer  where  2PA-generated  ESA  is  large  using  femtosecond 
excitation.  [54]  Figure  16  shows  excite-probe  data  for  two  excitation  energies  differing  by  a  factor 
of  2  along  with  fits  using  2PA  excitation  of  state  e  followed  by  ESA  to  state  u.  In  this  case,  there 
is  actually  a  multilevel  excited-state  system  where  the  2PA-excited  state  decays  to  a  lower-lying 
state  in  -240  fs  which  leads  to  the  delay  of  the  onset  of  ESA.  The  ESA  scales  as  the  input  energy 
squared  while  at  zero  delay  there  are  contributions  from  both  2PA,  which  scales  linearly  with  the 
energy,  and  ESA. 
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Figure  16.  Left  -  excite-probe  data  showing  the  change  in  optical  density  on  the  lead 
bis(ethynyl)porphyrin  molecule  shown  in  the  inset  right,  versus  temporal  delay  for 
excitation  energies  differing  by  a  factor  of  2.  Right  -  fits  for  the  two  energies  on  the  left 
using  a  combination  of  2PA  and  ESA.  Data  and  fits  provided  by  Joel  Hales  from  Ref.  [54], 

3.2.2  White-Light  Continuum  (WLC)  Excite-Probe  Spectroscopy: 

Excitation-probe  experiments  using  a  WLC  probe  as  shown  in  Fig.  17  have  greatly  speeded  data 
acquisition  of  nondegenerate  nonlinear  absorption  spectra  since  tuning  of  the  probe  becomes 
unnecessary.[55-57]  A  commercial  product  based  on  WLC  Excite-Probe  Spectroscopy  has  been 
developed  by  Helios  Inc. 

Usually  a  femtosecond  excitation  pulse  is  used,  which  is  also  used  to  create  a  femtosecond  WLC. 
Here,  the  probe  (WLC)  irradiance  is  kept  weak  to  avoid  any  nonlinear  effects  induced  by  the  probe 
itself.  Then,  in  most  cases  the  entire  WLC  can  be  incident  on  the  sample  at  the  same  time; 
however,  care  still  needs  to  be  applied  since  linear  absorption  of  the  WLC  can  produce  excited 
states.  Thus  the  femtosecond  excitation  prepares  the  sample  to  be  probed  by  the  WLC  which  can 
be  temporally  delayed  to  give  the  spectral  dynamics  of  the  NLA.  A  WLC  reference  can  also  be 
used  to  give  the  linear  transmittance  spectrum  which  can  be  used  to  determine  the  nonlinear  change 
in  transmittance  at  all  WLC  wavelengths  as  a  function  of  temporal  delay.[57] 
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Figure  17.  Femtosecond  excitation,  white-light  continuum  probe  spectroscopy  setup 
showing  a  2PA  measurement. [57] 

The  WLC  probe  may  be  used  to  measure  the  nondegenerate  2PA  spectrum,  a2(cop;  coe )  as  described 
by  Eq.  36.  Figure  1 8  shows  a  nondegenerate  nonlinear  absorption  spectrum  of  the  molecule  shown 
for  different  excitation  wavelengths  giving  different  levels  of  nondegeneracy.  For  the  largest 
nondegeneracy  the  nonlinear  response  is  enhanced  the  most  due  to  intermediate-state-resonance 
enhancement  (ISRE),  i.e.,  2PA  becomes  larger  if  the  intermediate  state  becomes  closer  to  being 
one-photon  absorbed. 
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Fig.  18.  Two-photon  absorption  spectrum  of  the  molecule  on  the  right  as  a  function  of  the 
photon  energy  sum  for  different  excitation  photon  energies  from  degenerate,  D  (blue 
circles)  to  very  nondegenerate,  ND  (green  squares). 


Additionally,  by  tuning  the  excitation  wavelength  to  a  linear  absorption  band  and  delaying  the 
WLC  probe  by  a  few  ps  with  respect  to  the  excitation  pulse,  we  can  measure  the  induced  absorption 
of  the  probe  over  a  broad  wavelength  range  to  obtain  the  ESA  spectrum.  Figure  19  shows  the 
ground  state  and  ESA  spectrum  of  a  thiacarbocyanine  dye  which  was  obtained  with  a  120  fs 
(FWHM)  excitation  pulse  at  900  nm  and  a  1 17  fs  (FWHM)  white-light  probe  pulse,  delayed  by  13 
ps  with  respect  to  the  excitation  pulse  over  the  spectral  range  of  400-800  nm.  In  this  case,  we 
observe  an  excited-state  cross-section  that  is  ~3  times  larger  than  that  of  the  ground  state  cross 
section  at  the  peak  spectral  position,  while  the  integrated  areas  of  their  absorption  bands  are  nearly 
the  same. 


X  (nm) 


Fig.  19.  Ground-state  absorption  (1)  and  ESA  spectrum  (2)  for  the  molecule  shown  at  the  bottom.  [58] 

Unfortunately  there  are  other  experimental  artifact  signals  that  appear  in  this  technique  that  require 
subtraction  to  obtain  the  signals  of  interest.  References  [55,  56]  describes  how  to  do  this.  The 
resultant  excite-probe  data  can  be  fit  using  a  single  equation  that  depends  on  the  WLC  chirp,  the 
group-velocity  dispersion  (GVD),  the  excitation  irradiance,  and  the  sample’s  nondegenerate 
nonlinear  absorption.  The  GVD  and  the  WLC  chirp  can  be  measured  using  the  same  apparatus  of 
Fig.  17  by  adding  appropriate  polarizers  as  in  the  optical  Kerr  effect  measurement.  [59]  This  excite- 
probe  spectroscopic  method  is  useful  for  measuring  both  2PA  and  ESA  and  can  distinguish 
between  these  nonlinearities  since  two-beam  experiments  show  the  temporal  dynamics. 

If  the  nondegenerate  nonlinear  absorption  spectrum  obtained  is  broad  enough  to  cover  all  the 
spectral  features,  Kramers-Kronig  relations  can  be  applied  to  calculate  the  nondegenerate 
nonlinear  refraction  dispersion.  [9,  60,  61]  In  addition  to  the  WLC  probe  for  measuring  NLA,  a 
narrower  spectrum  ultrashort  probe  pulse  can  be  spectrally  resolved  to  give  remarkably  sensitive 
measurements  of  NLR  in  addition  to  giving  the  NLA.  Here,  the  reported  sensitivity  to  a 
nonlinearly-induced  phase  change  is  ~10'7  waves. [62,  63]  This  method  was  first  used  by 
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Mysyrowicz  to  measure  the  nonlinear  refraction  of  air  and  temporally  resolve  the  various  nuclear 
components  coming  from  N2,  and  O2.  By  spectrally  resolving  the  transmitted  probe  pulse,  the 
cross-phase  modulation,  XPM,  can  be  seen  as  a  time-dependent  shift  of  the  frequency.  [63] 

3.2.3  Degenerate  Four-Wave  Mixing,  DFWM 

Another  excite-probe  measurement  is  degenerate  four-wave  mixing,  DFWM,  but  in  this  case  uses 
two  excitation  beams.  For  short-pulse  nonlinear  measurements  a  standard  geometry  is  shown  in 
Fig.  20  where  all  beams  can  be  temporally  delayed.  The  excitation  beams  are  usually  the  forward 
and  backward  excitation  beams  which  are  counter-propagating,  while  the  probe  comes  in  at  a  small 
angle  with  respect  to  the  forward  excitation.  The  analysis  usually  involves  discussions  of  the 
multiple  gratings  produced  which  scatter  the  various  beams  in  the  phase-conjugate  or  signal 
direction  which  is  backward  propagating  from  the  probe.  In  addition  to  the  usual  grating  terms 
there  is  a  so-called  two-photon  coherence  term  which  also  results  in  signal.  [64-66]  The  PC  signal 
using  this  geometry  in  the  low  depletion  (small  signal)  limit  is: 
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and  the  spot  sizes  and  pulsewidths  of  the  PC  signal  are:  wP£  -  Wj-  +  wb  +  wp  and 
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z PC  =  Tj  +Tjj+  rp  respectively.  The  important  result  to  note  is  that  the  nonlinear  susceptibility 

phase  is  not  given;  thus  both  real  and  imaginary  parts  contribute  to  the  PC  signal  and  they  cannot 
be  separated  without  determining  this  phase.  Various  homodyne  and  heterodyne  techniques  have 
been  developed  but  often  add  a  great  deal  of  complication. [67,  68]  The  geometry  discussed  next 
can  simplify  this  problem. 


Figure  20.  DFWM  setup  for  pulsed  measurements.  F,B  and  P  refer  to  the  forward 
excitation,  backward  excitation  and  probe  beams  respectively.  PC  is  the  phase  conjugate 
beam  or  signal  which  is  detected. 
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Another  standard  geometry  is  the  so-called  boxcar  geometry  shown  in  Fig.  21  where  three  input 
beams  (four  where  homodyne  detection  is  used  to  determine  the  phase)  are  arranged  to  give  a 
PC/signal  beam  which  forms  the  fourth  comer  of  the  “box”.  [69]  If  the  phase  of  the  signal  with 
respect  to  the  excitations  can  be  measured,  this  allows  separation  of  the  absorptive  and  refractive 
nonlinearities.  These  beams  can  be  generated  using  diffraction  gratings  (2D  gratings  can  produce 
all  the  needed  beams  from  a  single  input  and  essentially  eliminate  the  need  for  phase  stabilized 
beams). [69,  70] 


Figure  21.  Boxcar  geometry  for  DFWM.  The  dashed  line  is  the  phase  matched 
signal/PC.  For  homodyne  detection  a  weak  beam  is  also  incident  along  this  direction 
to  obtain  phase  information. 

The  magnitude  of  the  nonlinear  susceptibility  can  be  determined  with  the  use  of  a  known  reference 
sample  from; 
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where  L  is  the  sample  length,  subscripts  s  and  r  refer  to  the  sample  and  reference  respectively,  and 
T  is  the  linear  transmittance  of  the  samples.  Here  the  reference  sample  is  assumed  to  be 
nonabsorbing  and  the  experimental  conditions  identical.  The  signals  from  the  reference  and 
sample  are  IT  and  Is  which  are  the  maxima  in  the  DFMW  signals  with  time  delay.  Figure  22  shows 
an  example  of  a  DFWM  signal  in  the  boxcar  geometry  without  homodyne  detection  on  a  sample 
of  carbon  disulfide,  CS2,  where  the  solid  line  comes  from  a  calculation  of  the  expected  signal  from 
the  known  response  function  determined  in  Ref.  [70], 
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Fig.  22  DFWM  data  as  a  function  of  the  probe  time  delay  using  the  boxcars  geometry 
on  CS2  for  88  fs  and  42  fs  pulse  widths  of  the  excitation  and  probe  pulses,  respectively 
at  700nm.  The  solid  line  is  a  fit  using  the  bound  electronic  and  nuclear  responses  as 
discussed  in  Ref.  [70]  Inset  using  a  logarithmic  vertical  axis.  Taken  from  Ref.  [70] 

3.2.4  Two-Photon  Absorption  Induced  Fluorescence  Spectroscopy 


Many  organic  molecules  show  significant  fluorescence  quantum  yields.  This  allows 
measurements  of  both  one  photon  and  two-photon  absorption  induced  fluorescence  (2PF)  to  obtain 
a  degenerate  2PA  spectrum.  A  typical  experimental  setup  is  shown  in  Fig.  23. 


fluorophore 


Figure  23.  Experiment  schematic  for  the  two-photon  fluorescence  measurements 

After  simultaneous  absorption  of  two  photons,  the  excited  state,  rapidly  relaxes  to  the  zero 
vibrational  level  of  the  1st  excited  state).  The  fluorescence  intensity  is  a  measure  of  the  population 
of  the  first  excited  state,  which  is  proportional  to  the  2PA  coefficient.  Since  the  2PA  is 
proportional  to  the  square  of  the  input  irradiance,  /,  the  resulting  fluorescence  signal  is  also 
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proportional  to  I2.  This  dependence  is  thus  a  standard  criterion  to  distinguish  1PF  from  2PF  which 
was  first  used  by  Kaiser  and  Garrett  in  1961,  [71]  just  one  month  after  second  harmonic  generation 
was  reported.[72]  This  technology  was  further  developed  to  measure  2PA  spectra  [73,  74],  and  for 
two-photon  microscopy.  [75,  76]  Due  to  the  uncertainty  of  various  experimental  parameters,  e.g., 
collection  geometry  and  efficiency,  absolute  calibration  of  the  fluorescence  signal  to  determine  the 
2PA  is  problematic.  Standard  reference  molecules  with  known  2PA  spectra  are  often  used  to 
measure  relative  optical  fluorescence  signals. [73]  When  using  2PF  to  determine  2PA  spectra  as 
shown  in  Fig.  23  one  needs  to  minimize  strong  re-absorption  of  the  fluorescence  signal  for  high 
sample  solution  concentrations  by  focusing  the  excitation  beam  near  the  wall  of  the  cuvette  close 
to  the  collecting  lens.  The  fluorescence  spectrum  also  needs  to  be  corrected  for  the  spectral 
sensitivity  of  the  detector  and  for  any  residual  re-absorption  by  the  solution.  To  determine  the 
2PA  spectrum  of  the  sample  solution,  the  integrated  fluorescence  signal  is  then  compared  to  a 
reference  molecule  with  known  2PA  cross  section  under  identical  measuring  conditions.  The  2PA 
cross  section  for  the  sample  can  then  be  calculated  from: 
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where  <F>  is  the  integrated  fluorescence  and  <P>  is  the  average  excitation  power  which  is 
proportional  to  pulse  energy,  Cm  is  the  concentration  of  the  solution  measured  by  measuring  the 
1PA  spectrum  and  &f  is  the  fluorescence  quantum  yield  measured  by  1PF. 


There  are  several  other  factors  contributing  to  accurate  2PF  measurements.  The  concentration  of 
the  solution  should  be  as  high  as  possible  (but  below  the  aggregation  threshold)  to  increase  the 
probability  of  2PA,  thus  reabsorption  can  be  a  significant  issue.  Assuming  the  fluorescence 
spectrum  and  quantum  yield  are  independent  of  the  concentration  of  the  solution  and  the  type  of 
excitation  (i.e.  1PA  or  2PA),  the  2PF  spectrum  can  be  fit  to  the  1PF  spectrum  measured  at  lower 
concentration  (with  peak  OD  <  0.1)  to  obtain  a  correct  spectral  contour.  The  selection  of  a 
reference  molecule  is  based  on  the  spectral  range  of  the  fluorescence  of  the  sample.  A  large 
overlap  of  the  fluorescence  spectra  helps  minimize  the  measurement  error  induced  by  the  spectral 
sensitivity  of  the  detector.  Several  2PF  standards  have  been  published  with  emission  wavelengths 
ranging  from  390  to  1000  nm  and  excitation  wavelength  from  500  to  1600  nm.[73,  77,  78] 
However,  there  are  still  spectral  ranges  where  new/more  standards  are  needed. 


Figure  24  shows  the  2PF  spectrum  of  an  organic  dye  dissolved  in  chloroform  calibrated  against  a 
known  standard.  [78]  Z-scans  of  the  2PA  spectrum  are  also  shown  and  give  agreement  with  the 
2PF  data. 
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Figure  24.  Linear  absorption  spectrum  (black  line,  right  vertical  and  top  horizontal  scales), 
normalized  fluorescence  spectrum  (red  dash-dot  line,  top  horizontal  scale),  2PA  spectrum 
using  reference  calibrated  2PF  (green  symbols,  left  and  bottom  scales,  and  2PA  spectra  by 
Z-scan  (Black  squares)  of  the  molecule  shown  on  the  right.  From  Ref.  [78]. 

2PF  spectroscopy  is  in  principle  a  rapid  experimental  technique  to  measure  2PA  spectra  in 
comparison  to  other  commonly  adopted  techniques,  e.g.  Z-scan.  However,  verifying  the  quadratic 
dependence  on  irradiance  at  each  wavelength  considerably  lengthens  the  experimental  time.  It  is 
also  not  sensitive  to  the  ESA  that  can  follow  2PA,  since  the  ESA  doesn’t  contribute  to  the 
fluorescence.  The  sensitivity  of  2PF  can  be  very  high,  thanks  to  the  high  responsivity  of  PMT 
detectors  and  the  fact  that  2PF  is  a  zero  background  signal  technique.  A  2PA  cross  section  of  ~ 
1  GM  can  be  measured  given  a  sample  solution  with  large  fluorescence  quantum  yield  at  high 
concentration. 

However,  there  are  molecules  whose  fluorescent  quantum  yield  is  too  low  to  measure  by  2PF.  It 
is  also  difficult  to  measure  2PA  spectra  of  solid-state  samples  since  reference  material  geometries 
will  be  difficult  to  match  although  the  relative  spectral  shape  of  2PA  spectra  can  be  determined 
Finally,  the  2PF  technique  can  only  measure  2PA,  i.e.  imaginary  part  of  x(3>-  The  nonlinear 
refractive  index,  m,  has  to  be  measured  by  other  experimental  techniques. 

3.2.5  Fluorescence  Anisotropy 

One-photon  excitation  fluorescence  anisotropy  measurements,  especially  when  linked  to  quantum- 
chemical  calculations,  can  reveal  the  spectral  positions  and  orientations  of  the  transition  dipole 
moments  from  the  ground  state  So  to  the  first  excited  state  Si  (poi)  and  to  higher  excited  states  Sn 
(pon)  relative  to  the  orientation  of  the  emission  dipole  moment  (pio) .  This  information  cannot  be 
fully  obtained  from  linear  absorption  spectra.  The  excitation  anisotropy  spectrum,  r(Z),  is  defined 
as, 
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(41) 
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where  I\\  and  I|±  are  the  fluorescence  intensities  polarized  parallel  and  perpendicular  to  the 
excitation  light  respectively.  [79]  The  angle,  /?,  between  the  absorption  transition  moment  and  the 
emission  transition  moment  can  be  determined  from: [74] 
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Similarly,  2PA  anisotropy  can  be  measured  by  2-photon  excitation  and  measuring  the  ratio  of 
parallel  and  perpendicular  polarizations  which  for  a  3-level  system  gives  i*2pa  as: [80] 


18cos(y/2-/?01)cos/?01  cosy/2-7cos2  y/2  +  1 
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(43) 


Here  yis  the  angle  between  dipole  moments  participating  in  the  2PA  process,  and  j3o\  is  the  angle 
between  the  bisector  of  y  and  the  emission  transition  moment,  which  can  be  found  from  one- 
photon  anisotropy  measurements  in  Eq.  42. [81] 


Shown  in  Fig.  25  is  the  fluorescence  anisotropy  spectra  for  both  one  and  2PA  for  a  squaraine  dye 
along  with  other  data  described  in  the  next  section.  [82] 

Energy  (eV) 


Wavelength  (nm) 

Figure  25.  Shows  the  linear  absorption  spectrum  (blue  line)  of  the  squaraine  molecule  (SD 
2577  shown  in  Fig.  15)  with  the  1-photon  fluorescence,  1PF  (red  line),  1PF  anisotropy 
(green  circles),  2PF  anisotropy  (purple  squares),  and  the  2PF  and  Z-scan  measurements 
(red  squares)  of  the  2PA  spectrum.  [82]  Black  line  is  a  model  of  the  2PA  spectrum  based 
on  a  single  intermediate- state  and  three  2PA  states,  as  described  in  Ref.  [82], 
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4.  EXAMPLES  OF  USE  OF  MULTIPLE  TECHNIQUES 

In  this  section,  we  describe  how  organic  dyes  are  characterized  by  the  combined  use  of  several 
methods.  Figure  24  already  showed  a  dye  where  both  2PF  and  Z-scan  were  used  to  measure  the 
2PA  spectrum,  and  Fig.  25  shows  multiple  data  for  a  squaraine  dye  described  below. 

4.  1.  Squarine  Dye 

Figure  25  shows  the  results  of  both  linear  and  nonlinear  spectroscopy  of  a  Squaraine  dye  first 
reported  in  Ref.  [82],  The  linear  absorption  was  measured  in  solution  in  ethanol  using  a 
spectrophotometer  revealing  a  peak  absorption  wavelength  at  632  nm  with  a  peak  extinction 
coefficient  of  3.1  x  10-5  M_1cm_1.  The  fluorescence  spectrum  and  fluorescence  anisotropy 
were  measured  in  dilute  solution  using  a  spectrofluorimeter  from  PTI,  Inc.  As  described  above, 
the  one-photon  excited  fluorescence  anisotropy  reveals  the  position  of  the  one  and  two-photon 
transitions,  which  usually  cannot  be  determined  from  the  linear  absorption  spectrum. [79]  Also 
shown  in  Fig.  25  is  the  2-photon  excited  anisotropy,  which  shows  a  constant  value  at  all 
wavelengths  measured.  The  2PA  spectrum  shown  is  obtained  by  a  combination  of  2PF  and  Z- 
scan  measurements,  where  the  2-photon  fluorescence  is  measured  at  every  wavelength  and  in  each 
case  the  square-law  dependence  of  fluorescence  on  excitation  energy  was  verified  while  the  Z- 
scan  was  measured  at  three  wavelengths  to  provide  an  absolute  calibration  of  the  magnitude  of  the 
2PA  cross  section. 

The  experimentally  observed  2PA  spectrum  for  SD  2577  includes  3  bands:  a  weakly  allowed  band 
at  the  vibronic  shoulder  of  the  So— >Si  transition,  a  more  intense  band  at  A  =  408  nm  with  S2PA  = 
760  GM ,  and  a  much  more  intense  band  at  A  =  350  nm  with  82pa  =  5200  GM .  The  2PA 
spectrum  for  this  molecule  can  be  directly  mapped  onto  the  anisotropy  function  r(A) .  The 
observed  dips  in  anisotropy  at  410  nm  and  at  345  nm  correspond  to  the  peak  positions  of  the  2PA 
bands.  The  last  dip  at  309  nm  corresponds  to  a  position  near  twice  the  energy  of  the  main  So— >Si 
absorption  band.  As  a  result,  the  2PA  cross  section  at  this  energy  (corresponding  excitation 
wavelength  is  618  nm)  cannot  be  experimentally  determined  due  to  linear  absorption. 

4.2  Tetraone  Dye 

The  characterization  of  this  dye,  labeled  TD2765  shown  in  Fig.  26,  is  similar  to  the  squaraine  dye 
described  above.  The  methodologies  are  similar  to  those  described  in  Sec.  4.1  and  the  results  are 
shown  below  in  Fig.  26. 
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Figure  26.  Linear  absorption  spectrum  (blue  line)  of  the  tetraone  molecule  (TD2765) 
shown  on  the  left  along  with  the  1 -photon  fluorescence,  1PF  (red  line),  1PF  anisotropy 
(green  circles),  2PF  anisotropy  (purple  squares),  and  the  2PF  and  Z-scan  measurements 
(red  squares)  of  the  2PA  spectrum.  Solid  black  line  is  a  model  of  the  2PA  spectrum  based 
on  a  single  intermediate-state  and  three  2PA  states.  Structure  of  TD2765  shown  above 
plot.  From  Ref.  [83]. 

The  2PA  spectrum  in  this  case  was  also  taken  using  a  combination  of  2PF  and  Z-scan  methods. 
However,  close  to  the  linear  absorption  edge,  (A  «  700  nm)  the  Z-scan  gives  2PA  cross  sections 
as  large  as  35,000  GM.  Since  this  is  very  close  to  the  linear  absorption  region,  one  must  be 
cognizant  that  ESA  from  excited  states  produced  by  1PA  can  easily  be  mistaken  for  2PA,  as 
described  in  Sec.  2.  The  best  means  of  determining  the  mechanisms  for  the  nonlinear  absorption 
in  such  cases  is  to  perform  excite-probe  measurements,  since  the  2PA  has  an  ultrafast  response 
while  the  ESA  will  have  a  decay  time  equal  to  the  fluorescence  lifetime,  typically  hundreds  of 
picoseconds  to  nanoseconds.  Figure  27  shows  the  results  of  excite-probe  measurements  of  TD2765 
with  a  probe  wavelength  of  670  nm  and  three  different  excitation  wavelengths.  Exciting  at  710 
nm,  (curve  1),  no  ESA  is  observed  as  verified  by  the  full  recovery  of  the  probe  signal  within  the 
pulse  temporal  length.  At  700  nm,  (curve  2),  excite-probe  measurements  show  a  small  fraction  of 
ESA  (signal  lasting  after  the  pulse),  indicating  mostly  2PA  evidenced  by  a  recovery  of  94%  of  the 
normalized  transmittance.  Significant  ESA  is  observed  with  an  excitation  wavelength  of  690  nm 
in  (curve  3).  At  690  nm,  the  linear  absorption  cross-section  is  1.6  xlO'18  cm2,  almost  three  orders 
of  magnitude  smaller  than  the  peak  absorption  cross-section,  yet  there  is  significant  ESA.  If  the 
ESA  was  not  accounted  for  in  these  measurements,  it  would  have  appeared  that  we  measured  a 
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2PA  cross-section  in  excess  of  35,000  GM,  and  the  errors  would  probably  be  significantly  larger 
even  closer  to  the  absorption  edge.  Analysis  of  the  excite-probe  data  allows  us  to  determine  the 
relative  contribution  of  ESA  and  2PA  and  correct  the  Z-scan  data  to  give  the  2PA  spectrum  shown 
in  Fig  26. 


03 

O 

c 

03 

E 

< r ) 
c 
03 


"O 

<D 

N 

ra 

E 


1.00 
0.95 
0.90 
0.85 
0.80 
0.75 
0.70 
0.65 
0.60 

-1.0  -0.5  0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0 

Delay  Time  (ps) 


2 


■ 

t  * 

•  ■ 

▲ 

*1 

% 

1 - 1 - 1 - 1 - 1 - ' - 1 - 1 - 1 - 1 - 1 - 1 - 1 - » - 1 - 1 - 1 - 1 - 


Figure  27.  Normalized  transmittance  of  TD2765  as  a  function  of  probe  temporal  delay 
with  a  probe  wavelength  of  670  nm  and  3  different  input  excitation  wavelengths  getting 
closer  to  the  linear  absorption  edge  from  (1)  710  nm,  (2)  700  nm  to  (3)  690  nm  for  the 
tetraone  molecule.  The  absorption  peaks  at  625  nm  as  shown  in  Fig.  26.  This  shows  how 
resonant  excitation  of  excited  states  grows  to  nearly  dominate  the  absorption  as  the 
wavelength  is  reduced.  From  Ref.  [83]. 

To  measure  the  ESA  spectrum  shown  in  Fig.  28,  a  femtosecond  excitation  near  the  linear 
absorption  resonance  is  used  along  with  a  white  light  continuum  (WFC)  probe.  A  description  of 
this  technique  can  be  found  in  Refs.  [53,  58]  and  is  discussed  in  Sec.  3.2.2. 


Figure  28.  Absorption  cross-section  for  TD  2765  (3),  and  the  corresponding  excited-state 
absorption  spectrum  3’. 
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5.  OTHER  METHODS: 

In  this  chapter  we  have  discussed  many  methods  that  are  used  in  nonlinear  optical  materials 
characterization  along  with  other  selected  techniques.  There  are  several  other  methods  not 
discussed.  Here  we  mention  a  few  of  these  and  give  some  references  for  further  investigation. 
Thermal  lensing  from  2PA  [84-86] 

A.  Related  techniques  to  photothermal  lensing  [87,  88]  are: 

a.  Photothermal  displacement  [89] 

b.  Photothermal  beam  deflection  [90] 

c.  Photoacoustic  (optoacoustic)  detection  [91-94] 

d.  Thermal  calorimetry  of  2PA  [95] 

e.  Photothermal  interferometry  [96] 

B.  Beam  deflection.  [97]  This  is  a  new  excite-probe  method  that  that  allows  a  direct 
measurement  of  ultrafast  or  cumulative  nonlinear  refraction  and  nonlinear  absorption  to  be 
measured  with  a  demonstrated  sensitivity  to  nonlinear  phase  distortion  of  7720,000. [70] 

C.  two-photon-induced  phosphorescence  to  visualize  singlet  oxygen  concentration  [98] 

D.  Lifetime  imaging  to  differentiate  phosphorescence  from  fluorescence  [99] 

E.  Multiphoton  Ionization  Spectroscopy  [100-102] 

6.  CONCLUSION: 

There  are  many  different  experimental  methodologies  for  measuring  nonlinear  absorption.  The 
associated  nonlinear  refraction  can  often  lead  to  problems  in  unambiguously  determining  the 
nonlinear  loss;  thus  it  can  be  useful  to  simultaneously  determine  the  nonlinear  refraction.  This  can 
also  be  valuable  for  determining  the  underlying  physical  process  or  processes  which  can  often  be 
difficult.  Parametric  studies  of  the  wavelength  dependence,  temporal  dependence  and  sometimes 
the  optical  geometry  dependence  can  also  be  important.  The  variety  of  nonlinear  processes 
occurring  in  different  materials  is  rich  and  complicated,  and  it  is  rarely  the  case  that  a  single 
measurement  at  a  single  wavelength  with  a  single  pulsewidth  can  determine  the  physics  involved 
in  a  nonlinear  interaction.  Modem  nonlinear  spectroscopic  techniques  are  essential.  Ultrafast 
nonlinearities  often  require  femtosecond  resolution  to  unambiguously  determine  their  sign  and 
magnitude.  Longer  pulses  display  other  effects  such  as  excited-state  absorption  and  refraction 
which  can  even  dominate  the  nonlinear  response.  The  nonlinearities  occurring  in  organic  materials 
mimic  the  nonlinearities  in  semiconductors  where  free-carrier  effects  take  the  place  of  excited- 
state  effects.  However,  unlike  in  semiconductors  where  the  band  structure  determines  the 
response,  the  individual  molecular  levels  in  organics  can  be  engineered  allowing  the  nonlinear 


38 


response  to  be  altered.  However,  the  magnitudes  of  the  nonlinearities  are  still  not  as  large  as  is 
needed  for  many  of  the  applications  including  all-optical  switching  and  optical  limiting. 
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